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ABSTRACT 

A generalized wall function was proposed by 
Shih et al (1999). It accounts the effect of 
pressure gradients on the flow near the wall. 
Theory shows that the effect of pressure gra- 
dients on the flow in the inertial sublayer is 
very significant and the standard waU function 
should be replaced by a generalized wall func- 
tion. Since the theory is also valid for boundary 
layer flows toward separation, the generalized 
waU function may be applied to complex turbu- 
lent flows with acceleration, deceleration, sepa- 
ration and recirculation. This paper is to ver- 
ify the generalized wall function with numerical 
simulations for boundary layer flows with vari- 
ous adverse and favorable pressure gradients, in- 
cluding flows about to separate. Furthermore, a 
general procedure of implementation of the gen- 
eralized wall function for National Combustion 
Code (NCC) is described, it can be applied to 
both structured and unstructured CFD codes. 

INTRODUCTION 

Standard wall function (Millikan, 1938, Spald- 
ing, 1961, Tennekes and Lumley, 1972) has been 
widely and successfully used in various CFD 
applications for high Reynolds number turbu- 
lent flows. It is realized that the derivation 
of the standard wall function is valid only for 


flows with zero or small pressure gradients; how- 
ever, experienced CFD engineers can cdways 
find their ways to fix the problems where the 
standard wall function is not valid and does 
cause the troubles in numerical calculations. Of 
course, the degree of success varies among the 
individuals and the type of flows. Tennekes and 
Lumley (1972) described a method to deal with 
a boundary flow with zero wall stress. Shih 
et al (1999) extended this method to a general 
boundary flow and obtained a generalized wall 
function with the standard w'all function as its 
special case of zero-pressure gradient. The gen- 
eralized wall function is valid for various pres- 
sure gradients and zero-wall stress, hence over- 
comes the weakness of the standard waU func- 
tion. This will provide CFD engineers a unified 
way of using wall function for the CFD applica- 
tions without ad hoc fixing. We will first intro- 
duce the generalized wall function followed by 
the comparisons between the theory and numer- 
ical calculations of boundary layer flows with 
various pressure gradients. Then we describe 
the method of implementation of the wall func- 
tion for general 3-D CFD codes, especially for 
NCC code which is a 3-D unstructured code. A 
backward-facing step flow and dumpcombustor 
flows with and without combustion will be pre- 
sented as the examples. 
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GENERALIZED-WALL 

FUNCTION 


The generalized wall function in the inertial sub- 
layer proposed by Shih et al (1999) is listed as 
follows 
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( 1 ) 


where Ua is the mean velocity vector parallel 
to the wall or called tangential mean velocity, 
is the wall stress vector, (dP/dxi)tw is the 
streambned tangential pressure gradient at the 
wall. Ur and Up are the characteristic velocity 
defined by the wall stress and pressure gra- 
dient (dP/dxi)t^ as follows 


which shows that the mean velocity is no longer 
linear versus the normal wall distance due to 
the effect of pressure gradients. For the re- 
gion between the viscous sublayer and the in- 
ertial sublayer, i.e., the buffer layer, Shih et al 
derived a formulation following Spaldipg’s ap- 
proach. This formulation actually covers the 
whole surface layer from the wall all the way to 
the inertial sublayer. Therefore, they obtained 
an unified wall function valid for viscous sub- 
layer, buffer layer and inertial sublayer. See the 
paper of Shih et al (1999) for the details. 

BOUNDARY LAYER FLOWS 
WITH VARIOUS PRESSURE 
GRADIENTS 






u / dP\ 
p \ dx{ ) 


, 1/3 


Boundary layer flow with zero pressure 
^ ^ gradient 


which may be called the skin friction velocity 
and the pressure gradient velocity, respectively. 
Uc is a hybrid velocity scale consisting of Ur and 

Up . 


Uc = Ur + Up 


(3) 


This velocity scale Uc was used by Shih et al for 
scaling a general boundary layer flow. Note that 
Uc will never become zero by its definition, y is 
the normal distance to the wall, k is a constant 
0.41. Cl and C 2 are defined below 
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(5) 


where a ~ 5, /? « 8 and C « 5.0. Eq. (1) 
is valid in the region of inertial sublayer where 
>> 1 and y/6 << 1. For the region very 
close to the wall, u^y/u < 5 (viscous sublayer), 
Shih et al suggest 


Uit _ , 1 {dP/dx,)tu, f u^yY 

Uc pu"^\ 1/ J 2p uf \ P ) 

( 6 ) 


The calculations start with the uniform pro- 
files of mean velocity, turbulent kinetic energy 
and its dissipation rate. Our calculations with 
various grid resolutions, from 5 grid points to 
100 grid points, all approach the fuUy devel- 
oped turbulent boundary layer. Fig. 1 shows 
the skin friction versus the Reynolds number 
based on the momentum thickness. The ref- 
erence Reynolds number is 10000. Five calcu- 
lations with five different grid resolutions (uni- 
formly distributed across the boundary layer) 
are shown in Fig. 1 compared with the experi- 
mental data. Fig. 2 shows the results of mean 
velocity compared with the experimental data 
at Rgff = 7700. All calculations agree well with 
the experimental data when the calculated tur- 
bulent boundary is fully developed from their 
initial uniform profiles. Fig. 3 shows the cal- 
culations for different reference Reynolds num- 
bers, ranging from 10^ to 10^. The mean ve- 
locity profiles are plotted at Reg = 7700 and 
compared with experimental data. Good agree- 
ments are obtained. We conclude that the wall 
function implemented in this code has produced 
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Figure 1: Skin friction of boundary layer with 
zero pressure gradient 

the grid independent and reference Reynolds 
number (based on the arbitrary reference veloc- 
ity and length scale) independent solutions for 
fully developed turbulent boundary flows. 



Figure 2: Mean velocity of boundary layer with 
zero pressure gradient 

Boundary layer flow with adverse pres- 
sure gradients 

The wall function newly introduced by Shih et al 
(1999) includes the effect of pressure gradients. 
This waU function is claimed to be valid for large 
adverse pressure gradient including flows about 
to separate from the wall. To test this new wall 
function, we use a parabolic code and impose 
a constant adverse pressure gradient on a fully 
developed zero pressure gradient boundary layer 



Figure 3: Mean velocity of boundary layer with 
zero pressure gradient 

flow to obtain a fully developed boundary layer 
flow with the given adverse pressure gradient. 
Fig. 4 shows four cases, from the zero pressure 
gradient to the majcimum adverse pressure gra- 
dient. The latter corresponds to the boundary 
layer about to separate from the waU. The ve- 
locity scale in the figure Uc is defined as Ur + Up. 
Ur is the normal skin friction velocity, and Up 
is defined by the pressure gradient at the wall: 
Up = \{v j p)dP I dx\]l^ . Note that Uj and Up may 
individually approach to zero in certain cases, 
but their sum, i.e., u^ wiU never be zero in a 
boundary layer flow. The ratio, Up/uc, varies 
from zero to one. For zero pressure gradient, 
Up/uc = 0. For the flow about to separate from 
the wall, Ur must vanish and Up/uc = 1- Fig. 4 
shows the excellent agreement between the law 
of the waU (new wall function) and the compu- 
tations. Both the theory and the computations 
show the changes of the slope of the mean ve- 
locity profile as the adverse pressure gradient 
changes. 

Boundary layer flow with favorable 
pressure gradients 

For accelerating boundary layer flows, i.e. the 
flows with favorable pressure gradients, we use 
the same code and impose a constant favor- 
able pressure gradient on a fully developed zero 
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Figure 4: Mean velocity of boundary layer with 
adverse pressure gradients 



Figure 5: Mean velocity of boundary layer with 
favorable pressure gradients 


pressure gradient boundary layer to achieve a 
fully developed boundary layer flow with a given 
favorable pressure gradient. For an accelerat- 
ing flow, the skin friction is nonzero and in- 
creases with the pressure gradient, so that Up/uc 
is always less than one. Upjuc becomes small 
when Reynolds number is large. For example, 
in a fully developed channel flow (always fa- 
vorable pressure gradient), it can shown that 
y-phc = 1/(1 + where Re = Umaxh/u, 

Umax and h are the centerline mean velocity 
and the half width of the channel. Therefore, 
if Re > 1000, then Upjuc < 0.1. For accelerat- 
ing boundary layer flows Up/uc may be larger. 
Here, we have carried out the calculations from 
= 0 up to Upjuc = 0.3. The latter repre- 
sents an extremely large favorable pressure gra- 
dient under which the flow still remains turbu- 
lent. Five cases have been studied. The results 
are shown in Fig. 5. Again, the agreement be- 
tween the theory and computations is excellent. 
The mean velocity profiles will change signif- 
icantly with favorable pressure gradient if the 
Reynolds number is not very large. 

METHOD OF WALL FUNC- 
TION IMPLEMENTATION 

National Combustion Code (NCC) developed in 
NASA Glenn Research Center is a 3-D finite vol- 


ume unstructured CFD code. AH the variables 
are Cartesian. We will describe the method 
of implementation of the wall function for the 
NCC code. The method is general and can be 
extended to any other 3-D codes. In Cartesian 
coordinate system, the momentum equation can 
be written as 

dpUj dpu,Uj d _ dP 

dt dxj dxj ^ 

The integration of in an element (see 

Fig.6) gives 




kz=l 


k=l 


Av 


+ At; 

k=zl 


( 8 ) 


where Ay^k, -Az,k are the three compo- 

nents of face of the element. For the el- 
ements on the solid waJl, the values of nine 
stress components r,y at the wall (only six of 
them are independent) are needed, which can- 
not be accurately calculated using the values 
of velocity at the surrounding grid points with 
coarse numerical grids since the rapid change of 
turbulent mean velocity near the wall requires 
very fine grids to numerically calculate the wall 
stress. The wall function, however, can provide 
these wall stresses without numerical (differen- 
tial) calculations. 
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( 10 ) 



Wall 


Figure 6: Element on the wall 


Stress tensor in curvelinear coordi- 
nates 

Let us denote (x,y,z) as the Cartesian coordi- 
nates, 3.S the curvelinear coordinates. 

By definition and chain rule, we may write 


Tij = (m + Mr) 


d^n dx^ j 


( 9 ) 


. , (dr) dV 

r„ = 2(^ + MT)^-^ 


df] dVF\ 

dy dr] ) 


Eq.(lO) indicates that at the wall will be de- 
termined by the three derivatives of ^ and 
^ at the wall. The next section shows that 
these three derivatives can be provided by the 
turbulent wall function. 



Figure 7: Local curvelinear coordinates on the 
wall 


where Ui{L\V,W) are the velocity components 
in the Cartesian coordinates. Now, consider a 
wall surface element in the curvelinear coordi- 
nates. Its normal direction lies in 77, the surface 
is in the plane of ^ and ^ (see Fig. 7). At the 
wall, the derivatives of mean velocity with re- 
spect to the curvelinear coordinates are all zero 
except the normal derivatives: ^ and 

Therefore, the six independent at the wall 
can be written as 


Wall function 

The wall function gives relationship between the 
wall stress and the tangential velocity Ua in 
the inertial sublayer, that is Eq. (1). It can be 
written as 

U it ^ (^dP / dxi^tm 1 p ( 11 ) 

Uc p j Fi 


. .drjdP 

Txi = 2(/i + /xr)^-^ 

, ^fdridU 

rxv = (M + MT) [^-^ + 

, ,fdr)dU 
Txz - (m + Mt) + 

Tyy - 2(/i + Mr)^^ 


dx drj ) 
dr) 

dx dr) J 


where 

El 

F2 



+ Ci 
-f C2 


( 12 ) 


Eq. (11) provides us the wall stress or the nor- 
mal derivative of the mean velocity at the waU 
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dr]]^) without to carry out numerical 
differentiation. That is, if we write 


Wt = W- (20) 


/d[f dV dW\ 

^ \dT}'' drj’’ drj ) .^ 

then we may write 

'du\ _fUt 
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Calculation of dTjfdx, drjjdy^ drjjdz 


( 21 ) 


where, Ut, Vt and Wt are the three Cartesian 
components of tangential velocity in the inertial 
sublayer, and (vP)tw,x, (vP)tw,y and {vP)tw,z 
are the three Cartesian components of stream- 
lined tangential pressure gradient. The first nu- 
merical grid point is usually located in the iner- 
tial sublayer for CFD with coarse grids. How- 
ever, if the first grid point is close to the waU 
such that Ucv/v < 5, then the above relation 
must be replaced by the following relations; 



Ut- 
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[Vt-^ivPUyy\ 
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(18) 
(19) 


Calculation of tangential velocity and 
pressure gradient 

Let us denote the normal vector of the surface 
element by (5x,-, its length is n which is the nor- 
mal distance of the first grid point to the wall. 
We may calculate the three components of tan- 
gential mean velocity and streamlined tangen- 
tial pressure gradient as follows 


Ut = U-{U,6x,Y-^ 

V, = V- 


Let Xi be the location of the first grid point, 
(xf)i the location of the wall surface, 5,- the 
three components of the surface element with 
the area of S (5j + Sy + = 5^). Then the 

normal direction vector of the surface element 
will be Si/S. The three components of the nor- 
mal distance n (the first grid point away from 
the wall) can be calculated as 


6x = {x- a:/)y, 
6z = (z-zf)^ 


h = {y-yf)^, 
( 22 ) 


therefore, the normal distance n will be 


= 6x^ -f 6y^ -f 

^ (a; - xffSl + {y- yffSl + (z - z ffSl 

52 

(23) 

From Eq. (23), we obtain 


d„!dx = (* - 

n 

(24) 

n 

(25) 

dn/dz = 

n 

(26) 


6 


With all the above relations, we wiU be able to 
calculate the six at the wall using Eq. (10). 


COMPLEX-TURBULENT 

FLOWS 

Backward-facing step flow 

A backward-facing step flow of KKJ (Kim et al, 
1980) was used to verify the implementation of 
the wall function in the NCC code. The results 
using the new wall function are compared with 
the results from the old wall function built into 
the NCC code in Fig. (8). Apparently, the old 
wall function implementation method in NCC 
is not quite correct. Fig. (9) shows the effect of 
pressure gradient built into the new wall func- 
tion on a backward-facing step flow. The results 
show that the new wall function with pressure 
gradient effects gives better prediction in terms 
of the reattachment length (7.0 of step height 
due to the data of the experiment). The calcu- 
lation gives 6.8 with the effect of pressure gra- 
dient, versus 7.9 without the effect of pressure 
gradient. 


Dumpcombustor flow 

Another validation case is the dumpcombustor 
flow. Fig. (10) shows the flow with the combus- 
tion. The stream lines, the contours of pressure, 
NO and OH were plotted. Fig. (11) compares 
the same flow with and without combustion. 
Fig. (12) compares the calculations with the 
results from a low- Reynolds number turbulence 
model built into NCC, which does not need the 
wall function. This confirms that the results 
from the new wall function are quite reasonable. 


Conclusion 

The theory of generalized wall function has been 
verified by numerical simulations of fuUy devel- 



Figure 8: Backward- facing step flow 



Figure 9: Backward- facing step flow: effect of 
pressure- gradients 


7 








Figure 10: Dumpcombustor flow 1 


Figure 12: Dumpcombustor flow 3 




Cold flow (new wall function) 
(10000 iterations) 






Pressure 




oped turbulent boundary layer flows with vari- 
ous pressure gradients. It shows that the stan- 
dard wall function should be replaced by the 
generalized wall function for general boundary 
flows with large pressure gradients. A method 
of implementing wall function for 3-D complex 
flows has been introduced. The examples of 
complex 2-D flows show the success of imple- 
mentation of the wall function in NCC. Further 
validation against 3-D complex flows will follow. 
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